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Abstract 

The Fenna-Matthews-Olson (FMO) complex - a pigment protein complex involved 
in photosynthesis in green sulfur bacteria - is remarkably efficient in transferring ex¬ 
citation energy from light harvesting antenna molecules to a reaction center. Recent 
experimental and theoretical studies suggest that quantum coherence and entangle¬ 
ment may play a role in this excitation energy transfer (EFT). We examine whether 
bipartite quantum nonlocality, a property that expresses a stronger-than-entanglement 
form of correlation, exists between different pairs of chromophores in the FMO complex 
when modeling the EET by the hierarchically coupled equations of motion method. We 
compare the results for nonlocality with the amount of bipartite entanglement in the 
system. In particular, we analyze in what way these correlation properties are affected 
by different initial conditions. It is found that bipartite nonlocality only exists when the 
initial conditions are chosen in an unphysiological manner and probably is absent when 
considering the EET in the FMO complex in its natural habitat. It is also seen that 
nonlocality and entanglement behave quite differently in this system. In particular, for 
localized initial states, nonlocality only exists on a very short time scale and then drops 
to zero in an abrupt manner. As already known from previous studies, quantum entan¬ 
glement between chromophore pairs on the other hand is oscillating and exponentially 
decaying and follow thereby a pattern more similar to the chromophore population 
dynamics. The abrupt disappearance of nonlocality in the presence of nonvanishing en¬ 
tanglement is a phenomenon we call nonlocality sudden death] a striking manifestation 
of the difference between these two types of correlations in quantum systems. 
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INTRODUCTION 


In recent years, there has been an increasing interest in studying biological systems in terms 
of the existence of nontrivial quantum effects®^. Especially the excitation energy transfer 
(EET) between chromophores in photosynthetic complexes has been heavily studied ever 
since quantum coherence between excited states in the Fenna-Matthews-Olson (FMO) com¬ 
plex was experimentally verihed by two dimensional electronic spectroscopy*^. Not only was 
the existence of coherent pathways in the FMO complex revealed, quantum coherence was 
also shown to last much longer than expected for such large and noisy systems, and a hy¬ 
pothesis is that this may play a role in the known very efficient EET in photosynthetic 
complexes*^. Since then, several attempts to explain why photosynthetic complexes could 
beneht from coherent EET have been made*®. One such is that the quantum coherence could 
help to direct the energy flow towards the reaction center in a unidirectional manner*®. 

The experimental evidence for coherent EET in the FMO complex motivates the develop¬ 
ment of new methods to model the excited state dynamics of chromophores in a surrounding 
protein scaffold. It has become clear that Markovian models cannot fully capture the be¬ 
haviour of the system-environment interactions in these systems and hence non-Markovian 
models have to be considered. Based on the hierarchial expansion technique proposed in 
Refs IMill, Ishizaki and Fleming rehned the theoretical framework and developed a tool for 
investigating the excited state dynamics in a photosynthetic complex*®^. This form of hi- 
erarchially coupled equations of motions (HEOM) has served as the benchmarking method 
for these systems since it is able to describe quantum coherent motion and incoherent elec¬ 
tron hopping in a unihed manner*!^. Using this method, Sarovar et al^ examined quantum 
entanglement in the FMO complex and found that bipartite entanglement of chromophores 
exists on a time scale relevant for the mechanism of EET. Since entanglement can be used as 
a resource in quantum information processing tasks, it was speculated if the hndings could 
play a role in explaining the near unity efficiency for these complexes to convert solar light 
into chemical energy. 

Quantum nonlocality, as introduced by Bell*^, is another correlation property of com¬ 
posite quantum systems. While quantum entanglement, already mentioned by Schrodinger 
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as the key property that distinguishes quantum mechanics from the classical world^, is 
a meaningful property only for quantum systems, nonlocality is given by a criterion that 
is formed without specifying whether the systems are quantum mechanical or not. In the 
quantum-mechanical context, although entanglement and nonlocality are the same for pure 
states, this is not in general true for mixed states, i.e., they are inequivalent properties for 
composite quantum systems^^^^. In particular, nonlocality is a stronger correlation property 
than entanglement in the sense that nonlocality implies entanglement but not vice versa^. 

The aim of this study is to examine whether bipartite nonlocality can exist during EET 
in the FMO complex, which would imply a striking nonclassical feature of this process. 
By studying the existence of quantum nonlocality in such a complex, it is possible to add 
insights into if and how quantum effects play a role in photosynthesis; insights that could be 
of interest for artificial photosynthesis and solar cells as well as for quantum computation. In 
recent studies^^^^, a considerable enhancement in the efficiency due to quantum coherence 
was found in model systems of photosynthesis mimicking solar cells. If such discoveries also 
can be connected to the existence of entanglement and nonlocality, it would be desirable to 
explore the underlying physics of such phenomena further. 

MODELING THE EET IN THE FMO COMPLEX 

EET in the FMO complex has been modeled by employing the HEOM method*^^. A 
description of the FMO complex as well as the quantities and conditions used to model 
EET, is given followed by a brief review of the HEOM based model of the FMO complex, 
modified by including an explicit mechanism of the excitation energy trapping at the reaction 
center. 


The FMO complex 

The FMO complex is a photosynthetic complex found in green sulfur bacteria, Chlorohacu- 
lum tepidum. It consists of three identical monomers, each containing seven chromophores 
embedded in a protein scaffold, as can be seen in Fig. The structure of this complex as 
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Figure 1: The FMO complex trimer and its position between the antenna molecules and the 
reaction center to the left, and the chromophores in their protein scaffold for one monomer 
to the right. The chromophores are chlorophylls in this particular system. 

well as its site excitation energies and inter-site couplings have been experimentally investi¬ 
gated by different spectroscopic methods^^^^^® and the system has been used as a model for 
larger light harvesting complexes for more than 35 year^^H 

The FMO complex does not include any antenna molecules. In other words, it is not 
responsible for capturing the light energy, but only for transferring it to the reaction center. 
The antenna molecules are located in such a way that they can transfer excitation energy 
to the three monomer units. We restrict the system studied to one monomer; a reasonably 
simplifying assumption in the description of correlations accompanying EFT in this system. 

The individual chromophores in the FMO complex can be modeled as two level systems 
(TLSs) by taking into account only the So —)• Si transition, with Sq and Si being the singlet 
ground state and singlet first excited state, respectively. Furthermore, since Green sulfur 
bacteria recieve very little sunlight in their natural habitat^, it is reasonable to assume that 
the FMO complex only contains one such excitation at a time. This reduces the electronic 
Hamiltonian of the FMO complex (one monomer) to have the form 

7 7 

k=l k^l 

where Ek is the electronic excitation energy of chromophore k when being in its excited 
state, while the other chromophores remain in their ground state. This corresponds to the 
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localized state 


\k) = \Sl..Sl..Sl), (2) 

where the superscript denotes the chromophore site. Furthermore, J^i is the dipole coupling 
describing the electrostatic interaction between the charge distribution of chromophores k 
and 1. It depends strongly on the dipole moment orientations and relative positions of the 
chromophores in the protein complex structure^. 

Calculations suggest that the most favourable structure is that chromophores 3 and 4 
are linked directly to the reaction center^. This result has been confirmed experimentally 
when examining the FMO complex with masspectrometry®. Since the structure of the FMO 
is known^, this implies that chromophores 1 and 6 connect to the base plate (that is the 
part of the antenna molecule complex that connects it to the FMO). The simplest way to 
model EFT in the FMO complex would hence be to assume a localized excitation on either 
chromophore 1 or 6 as the initial conditions for EET. This can be written as 

^focalized = I, a; =1,6. (3) 

However, since the distance between the antenna molecules and the base plate as well as 
the distance between the base plate and the FMO are so large (~2 nm^^MHI ~1.5 nm^^, 
respectively) in comparison to the intracomplex distances, a model that better captures the 
condition for initial excitation of the FMO complex in its natural habitat is to assume that 
the excitation is transferred from the base plate to the FMO by Forster resonance energy 
transfer (FRET)^. This would populate the FMO exciton states le^), being the eigenstates 
of the Hamiltonian in Eq. ([^, in proportion to their occurrence at chromophore 1 or 6. 
These FRET initial states can hence be written as 

7 

^FRET = X] l(2^lo)|^ |er-)(er.|, X = 1, 6. (4) 

r=l 

In this work, we shall examine correlations between chromophore pairs in EET arising from 
the initial states given by Eqs. ([^ and (|^. 
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The HEOM method for the FMO complex 


is a numerically exact method in which the environmental influence on a quantum 
system is treated in a statistical manner. For the FMO complex, the system of interest is 
the seven chromophores in one of the monomers with the Hamiltonian given in Fq. Q. The 
protein scaffold is modeled as a set of harmonic oscillator modes, i.e., a phonon bath, which 
represents nuclear motion, both intramolecular and those of the environment. The transfer 
of excitation energy from one chromophore to another occurs via nonequilibrium nuclear 
configuration (phonon states) in accordance with a vertical Franck-Condon transition. The 
phonons, locally and linearly coupled to each chromophore fc, relax to their equilibrium states 
under energy loss. This so-called reorganization energy, denoted Xk, can be measured via 
the Stokes shift^^EZl^ 

The memory of the local environment of chromophore k is characterized by a relaxation 
function Tk{t), modeled as an overdamped Brownian oscillator, which takes the form 

Tk{t) = 2Afce-^'=h (5) 


The parameter jk represents the time scale of the fluctuations and energy dissipation for 
the kth chromophore and is related to the non-Markovian behaviour of the dynamics. The 
relaxation function, and hence 7 ^, can be investigated by two dimensional electronic spec¬ 
troscopy*^^. 

The time dependent density operator g{t) describing FFT in the FMO complex is ob¬ 
tained by solving a set of hierchically coupled equations of motion, 

7 


d . 




y^^nk'^k 

k=l 


C(n,«) + ^ 4sC(ni!+. *) + t) 


( 6 ) 


k=l 


where the operator C{0,t) is identical to g{t) while the higher order operators C(ii 7 ^ 0 )^) 
are auxiliary operators. Here, n is a set of nonnegative integers, n = (ni, n 2 ,.., n 7 ). The 
notation nfc+(nfc_) stands for adding (substracting) 1 to the corresponding Uk in n. The 
Liouvillian superoperator Cg is composed of the electronic Hamiltonian of the FMO complex 

and the site reorganization energies, and takes the form 

7 

Cgg = [He + y2^k\k){k\,g] (7) 

k=l 


6 






with [•, •] the commutator and g any linear operator acting nontrivially on the single exci¬ 
tation subspace of the full Hilbert space of the FMO complex. The superoperators and 
are so called phonon-induced relaxation operators and correspond to the influence of the 
environmental fluctuations. They have the form 

^kg = i[Vk,g] (8) 

and 

Qkg = i‘^[Vk,g] +^7k{Vk,g} (9) 

with Vk = \k){k\ and {•, •} the anti-commutator. Here, f3 = {kBT)~^ with ks being the 
Boltzmann constant and T the temperature of the phonon bath. 

The HEOM can be truncated by setting 

^C(n,t) =-*>CeC(n,t), (10) 

for all n satisfying ni + n 2 + ... + n-j = N, N being the truncation level. This condition will 

terminate the generation of auxiliary operators. 

Following Ref.l^, a Liouvillian that models the trapping of excitation of chromophores 3 
and 4 is added to the right-hand side of Eq. ([^. It has the form 

ArapC(n,t) =-rtrap{|3)(3|,C(n,f)}-rtrap{|4)(4|,C(n,f)}, (11) 

where rtrap is the trapping rate, assumed to be the same for chromophores 3 and 4. When 
comparing to Eq. (§, it can be seen that adding this Liouvillian makes the population of 
chromophore 3 and 4 decay faster than for the other chromophores. 

The density operator describing EET in the FMO complex can be written as 

7 7 7 

Q{t) = + + 4i(^)IO(^l), (12) 

k=l k=l l>k 

where \k) are the site basis states defined in Eq. (|^. Here, Qkk is the population of an 
excitation at chrompohore k and Qki describes the coherence between chromophores k and 1. 
Note that g can be viewed as a 7 x 7 dimensional Hermitian matrix due to the restriction 
to one coherent excitation at each instant of time. 
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MEASURE OF BIPARTITE QUANTUM CORRELATIONS 


Quantum nonlocality is a property of composite quantum systems whose subsystems show 
correlations that are too strong to be explained by a local realistic theorj^, i.e., a theory 
where physical variables (such as, e.g., positions and momenta of particles) are assumed to 
have well-defined local values prior to measurement. The existence of quantum nonlocal¬ 
ity was discovered when Bell derived^ an upper bound (Bell’s inequality) for correlations 
between two systems to be local, and then showed that the correlations within certain com¬ 
posite quantum systems may exceed this limit. 

Since Bell’s original formulation of his inequality, there have been other Bell-like in¬ 
equalities. One such that is suitable for investigating bipartite nonlocality for TLSs like the 
chromophores in the FMO complex is the Clauser-Horne-Shimony-Holt (CHSH) inequality®. 
It states that any set of variables a, a', b, and b' that can take values ±1 must satisfy 

\{ab) -|- {ab') + {a'b) — {a!b') \ < 2 (13) 


if a local realistic theory applies to the pairs a, a' and b, b' at distant locations. Here, {ab) 
denotes the average of the product of the outcomes a and b. 

To test the validity of the CHSH inequality for measurements on two distant TLSs, each 
characterized by the Pauli operators a = (di, <72, da), described as a composite quantum 
system with density operator p, the operator 

®CHSH = (a • d) (g) (b ■ d) -h (a ■ d) o (b' ■ d) 

-h(a' • d) g (b • d) — (a' ■ d) g (b' ■ d) (14) 


with a, a', b, and b' unit vectors in can be used, as the measurements of (a-d), (a'-d), (b- 
d), and (b'-d) on the respective TLSs have outcomes ±1. Thus, by comparing with Eq. ( |T^ , 
one may conclude that the correlation in p is nonlocal, i.e., cannot be accounted for by any 
local realistic theory, if there exists a, a', b, and b' such that | (IBchsh) I = |Tr ^pfScusH j I > 2. 

A necessary and sufficient condition for the correlation between any two TLSs to be 
nonlocal has been developed by Horodecki et a/.®. This criterion is based on a quantity 


M(/3) that maximizes the expectation value of the Bell operator in Eq. (14) such that 


max |(iBcHSH)| = 2^M(p). 

a,a',b,b' 


( 15 ) 




M(p) is found as 


■^(p) Pp l-^qi 


(16) 


where /ip and fiq are the two greatest eigenvalues of |T(/i)|^, T(p) being the correlation matrix 


T(p) = 


(17) 


til ti 2 tl3 
^21 ^22 ^23 
^^31 ^32 ^33 y 

with matrix elements = Tr(/3(Ja ® d/?) for all combinations of Pauli operators da- The 
correlation is nonlocal whenever 1 < M(/3) < 2 , where the maximum value is given by the 
Cirel’son bound^ K®chsh)| < 2 a/2 . This motivates that the quantit}^ 

B{p) = Amax{(M(p) - 1),0} (18) 

can be used as a measure of the amount of nonlocality. This Bell-CHSH measure is directly 
comparable with concurrence, defined a^^ 

C(p) = max{Ai - A 2 - A 3 - A4, 0}, (19) 

where Ai,... A 4 are the decreasingly ordered eigenvalues of the positive operator 

R = \J vA <^2 ® 


(J2P"(J2 (» CT2 


( 8 ) cr 2 vA 


( 20 ) 


with complex conjugation taken in the computational product basis. Concurrence uniquely 
determines the entanglement of format ion^S^l of two TLSs; as such, C{p) is an entanglement 
measure in its own right. It has been used to analyze the amount of bipartite entanglement 
in the FMO complex by Sarovar et al^. 

The time evolution of the composite system of chromophores n and m in the FMO 
complex can be found by calculating the reduced density operator from the full density 


operator q given in Eq. (12). This is done by tracing q over the other five chromophore 


degrees of freedom. The resulting reduced density operator takes the form 


k,l,p,q=0,l 

= [TT{g)-igmm + gnn)]\S^S^){S^S^\ 

I I cm C'R \ / cm qti I i ^ I cm on \ / cm on I 

' ^mm\^0 / Wo *^1 I t/mn Wo /\‘^1 *^0 I 

I * I om Qn\ j nm on I i ^ I cm cn\ j nm on \ 

“Tt/mnWl ^0/\^0 I t/nnWl ^0/\^l ^0 h 


( 21 ) 
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where we have taken into account that HEOM and excitation trapping are not trace pre¬ 
serving, i.e., that Tr(p) < 1. Note that pmn can be viewed as a 4 x 4 Herniitian matrix. One 
finds 

h"! I~^2 4 |Pmn| , ^3 [Tr(p) 4“ ^^nn)] ) 

which implies 

B{pran) = l^maX { (max{8 \Qmn\^ A\Qmn\^ + [Tr(p) - 2{Qram + Qnn)f} “ l) ,0}. (23) 

Similarly, concurrence for this class of states takes the form 

= 2|p^„|, (24) 

which coincides with the square root of pi. 

Recent attempts to quantify the amount of coherence in quantum states^^^^ have lead to 
different types of coherence measures. One of these is the li norm measure of coherence^, 
defined as 

Ci.(p) = 5^|ph|, (25) 

k^l 

i.e., the sum of the absolute value of the off-diagonal elements of the density operator. For 
Pmn this reduces to 

C'p^Pmn) 2 I Pmn I ; (^6) 

which implies that concurrence and the li norm measure of coherence are in fact identi¬ 
cal quantities for all chromophore pairs. Thus, a nonzero concurrence can equally well be 
interpreted as an expression for having a nonzero coherence, rather than being a sign of 
correlation. This is a consequence of the restriction to the single excitation subspace of the 
Hamiltonian in Eq. ([^. Thus, the restriction itself implies the existence of entanglement as 
a consequence of the nonvanishing coherence in the system. 

Nonlocality, on the other hand, is essentially different from coherence as it also depends 
on the populations Qmm and Qnn of the chromophores via the nondegenerate eigenvalue pa 
of |r(p)p. This makes the Bell-CHSH measure a proper quantifier of correlation between 
chromophore pairs. 
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We may ask under what circumstances a chromophore pair exhibiting nonvanishing co¬ 
herence is nonlocally correlated. As follows from the explicit expressions in Eq. (23), the 


nonlocality condition relies on the relation between the populations Qmm and Qnn of the pair 
and the coherence 2 Since the latter is bounded by the former as^^ 


I Qmn I ^ \/ Qmm Qnn j 


(27) 


which follows from positivity of the reduced density matrix, we see that /xi and /i 2 become 
very small unless the system is considerably localized to the (m, n) pair, in case of which 
/i 3 also becomes large. Alternatively, fis can be close to the critical value 1 sufficient for 
nonlocal correlations if Qrnm and Qnn are both very small (in the order of a few %, say) and 
Tr(^) < 1; however, this cannot give rise to any nonlocality. To see this, we hrst note that 


Af (Pmn) P-l -|- P-s ^ ^QmmQnn T (Tr(^) ‘^Qmm ‘^Qnn) 

= Tr(^) (Tr(^) - Agmm - 4^nn) + + 12^, 

^ 1 4 {^Qmm. T Qnn Qmm Qnn ^QmmQnn') 


mm Qnn 


(28) 


by combining Eqs. (22), (27), and Tr(^) < 1. It is straightforward to see that Qmm + Qnn 


q\^ — SQmmQnn > 0 for Small Qmm and Qnn, which implies that 


MiPmn) < 1 (29) 

excluding the possibility of nonlocal correlations in this case. We conclude that only chromphore 
pairs for which the population is large can be nonlocally correlated. 

NUMERICAL DETAILS 

The parameters of our numerical model are chosen in accordance with previous work on 
EET in the EMO complex. As the electronic EMO Hamiltonian for Chlorobaculum tepidum 
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in the chromophore site basis, we use^ 
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(30) 


The reorganization 

energy and the relaxation time-scale 7 *, are assumed to have the same values, 35 cm“^ and 
50 fs“^^, respectively, for all seven chromophore^. The time scale for the trapping by 
the reaction center is set to 1 p^lMI ^nd the bath temperature to 300 K (same as in Ref.^. 
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Figure 2 ; Convergence of the density operator with respect to the truncation level of the 
HEOM. The change in the density operator when going from level N to level A^-|-1 is charac¬ 
terized by the logarithm of the corresponding trace distance D(iV, N + 1) = D^g^, Qn+i) = 
|Tr Ipat — Pat+iI- The initial excitation is either on chromophore 1 (upper panel) or chro¬ 
mophore 6 (lower panel). Notice how at the higher truncation levels {N > 10) the trace 
distance is comparable to the numerical noise. 
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The HEOM given by Eq. ([^ is numerically integrated in the range 0 to 1000 fs using 
the Runge-Kutta-Dormand-Prince method® To measure the convergence of the HEOM 
solution, we use the trace distance® 

D{gN, Qn+i) = Iqn — ^w+il (31) 

for density operators and gN+i at truncation level N and iV + 1, respectively. In our 
simulations in the next section, we use a truncation level iV = 12, which implies an accuracy 
in the order of 10“®, as can be seen from Fig. 

RESULTS AND DISCUSSION 

Our simulations show that no nonlocality is found when the initial states are given by Eq. Q, 
corresponding to FRET from the antenna molecules to the FMO, although entanglement 
still exists. In particular, a considerable amount of quantum entanglement is found for 
chromophore 1 and 2 when the exciton states are projected on chromophore 1 as well as for 
chromophore 5 and 6 when the exciton states are projected on chromophore 6. These results 
can be seen in Figs. [^and|^ On the other hand, bipartite quantum nonlocality exists when 
the localized initial conditions according to Eq. ^ are used. For these initial conditions, 
nonlocality is found for two chromophore pairs; chromophore 1 and 2 when chromophore 1 
is initially excited, and chromophore 5 and 6 when chromophore 6 is initially excited. These 
results are presented in Figs. and together with the time dependence of the bipartite 
entanglement for the same pairs of chromophores. 

To get father insights into these Endings we analyze in more detail the structure of the 
evolution arising from the two types of initial conditions. Let us start with the FRET case, 
where the initial states are given by Eq. (j^. By writing the exciton states \er) = X]I=i Crfcl^), 
we find 

7 

2fRET = C^xkr)(er|, (32) 

r=l 

where we have used that all Crk are real-valued and we have ordered \er) with increasing 
energy. The reduced density matrix for the (m, n) pair is characterized by the site probablities 
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Figure 3: Nonlocality B{pi 2 ) and entanglement C{pi 2 ) for the chromophore pair 1 and 2 as 
a function of time for FRET initial state when exciton states projected on chromophore 1 is 
used as initial condition. Note that the nonlocality B{pi 2 ) vanishes at all times. 


Qmrrf, Qnni and the off-diagonal term g^n given by 

7 


E 2 2 

^rx^rmi 

r=l 

7 

E 2 2 

^rx^rn-i 

r=l 

7 

E 2 

r=l 


(33) 


We focus on the x = 1 case. As demonstrated above, only chromphore pairs for which 
the population is large can exhibit nonlocal correlations; thus, we only need to consider the 
dominant terms in the density operator of the FMO complex. By inspection of the explicit 
eigenvectors, we hnd that Pfret dominated by the pure state components jea) and jee). 
In fact, ~ 0.769 and Cg^ ~ 0.208; thus, these two exciton states populate almost 98% of 
this initial state and we can safely ignore all the other exciton states. We further find that 
both lea) and \e7) are to a large extent localized to the first two chromophores, as can be 
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Figure 4: Nonlocality B{p^q) and entanglement C{p^q) for the chromophore pair 5 and 6 as 
a function of time for FRET initial state when exciton states projected on chromophore 6 is 
used as initial condition. Note that the nonlocality B{pi 2 ) vanishes at all times. 

seen from the expansion coefficients 

C31 = 0.877, C32 = 0.440, 

C61 = -0.456, C62 = 0.871, (34) 

obtained by numerical diagonalization of the electronic Hamiltonian He. Thus, correlation 
is essentially concentrated to the first two chromophores, and can be expressed in terms of 
the matrix elements 


Q 12 ~ C31C32 + C61C62 

= IC 31 C 32 I — Cgg |C 6 lC 62 | = 0.214, 

Qii K. C 31 + Cgi = 0.635, 

Q 22 ~ ^ 3^032 + Cg 3 Cg 2 = 0.307. (35) 
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Figure 5: Nonlocality B{pi 2 ) and entanglement C{pi 2 ) for the chromophore pair 1 and 2 as 
a function of time when chromophore 1 is initially excited. Note how the nonlocality drops 
to zero within about 80 fs. 

These numerical values imply 

C{pi2) = 2|pi2| ^ 0.428, 

hi = h2 = [C'M'^ 0.183, 

h3 = (l-2pn-2p22)'~ 0.781; (36) 

thus, Pfret shows entanglement {C > 0) but no nonlocality (5 = 0 since Pp + Pg = 
0.183 + 0.781 = 0.964 < 1) between the hrst two chromophores. Since correlations essentially 
are concentrated to this chromophore-pair, it follows that all chromophore pairs are locally 
correlated at f = 0, given the FRET initial condition. 

On the other hand, we note that the two dominating exciton states p®® = |e 3 )(e 3 | and 
p®® = lee)(eel are strongly nonlocal between chromophores 1 and 2; indeed, one hnds 

Bipll) = C'(pf 2 )^ 2 |c 3 lC 32 |= 0 . 801 , 

Bipll) = W2)^2|ceice2|= 0.822, (37) 

where we have used that B = C for pure states of any TLS-paii®. By comparing these 
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Figure 6 : Nonlocality B{p^q) and entanglement C{p^q) for the chromophore pair 5 and 6 as 
a function of time when chromophore 6 is initially excited. Note how the nonlocality drops 
to zero within about 40 fs. 


expressions with the expression for qi 2 in Eq. (35), we find 


Q12 ~ ^ZlB{Pi2) 


(38) 


which explicitly entails that the essential source for the disappearance of nonlocality when 
mixing the exicton state is the destructive quantum-mechanical interference (relative minus 
sign) between the first-site components C 31 and cei of jea) and jee). 

A similar analysis of the x = 6 case can be carried out with the same conclusion that all 
chromophore pairs are locally correlated, but entangled at t = 0 . 

FRET are stationary states under the action of the electronic Hamiltonian, but they may 
undergo nontrivial evolution under influence of the environment. Thus, the environment is 
a potential source of nonlocal correlations to show up at t > 0. However, as can be seen in 
Figs. and it turns out that nonlocal correlations never appear given the FRET initial 
conditions, despite the fact that entanglement is typically nonvanishing. 

Next, we consider the case of localized initial states p^jcaiized = ^ = 1)6. These 

states are not eigenstates of Hg and can therefore in principle give rise to nonlocal correlations 
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during time evolution. 

Let us determine for which chromophore pair(s) nonlocality can occur. In the short time 
limit, we may expect that n = 0 is the dominant term in the HEOM. Thus, we may assume 
that g{St) ~ small 6t (from now on, we put h = 1 for notational 

simplicity). By expanding in 6t, we hnd to leading orders: 


C{p^n) ~ 26t\J, 


C(^Pmn) ~ 25t \JmxJxn \ 5 


B (^Pxn) 
B{Pmn} 



max 


0 , 



(39) 


where m,n ^ x. Here, Jki are the real-valued dipole coupling parameters given by the off- 
diagonal terms of Hf.. We may therefore draw the following conclusions about the short-time 
behavior of a localized initial state: 


(i) entanglement between the initially excited chromophore x and any other chromophore 
n increases linearly with 6t with proportionality factor given by the corresponding 
dipole coupling parameter \Jxn\', 

(ii) entanglement between any other pair of chromophores increases quadratically with 6t-, 


(iii) nonlocal correlation between chromophore x and chromophore n increases linearly with 
5t provided 


B > V 72. 

'^xn ^ ^ "^xli 

l^n 

(iv) nonlocal correlations between any other pair of chromophores vanish. 


(40) 


The inequality in Eq. (40) can be satished for at most one pair. The explicit form of 


Hg entails that only chromophores 1 and 2 can show nonlocal correlation after excitation 
of chromophore 1, while only chromophores 5 and 6 can show nonlocal correlation after 
excitation of chromophore 6 . For these pairs we further notice that entanglement grows 
linearly with almost the same speed since | J 12 I ~ As can be seen in Figs. and 
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these results are confirmed in our numerical simulations of the HEOM, where we find 
nonlocality for t less than ~ 80 fs and ~ 50 fs for x = 1 and x = 6, respectively. Indeed, 
at these instances of time, nonlocality suddenly disappears, although entanglement persists. 
We call this phenomenon nonlocality sudden deaths being a striking manifestation of that 
nonlocality and entanglement are different properties for mixed quantum states^^. 

The absence of nonlocal correlation after its sudden death for localized initial states 
is a consequence of open system effects. In contrast, the complete absence of nonlocality 
in the FRET case is a consequence of the local nature of the initial state. Thus, we may 
conclude that the disappearance of nonlocality has different origin for the two types of initial 
conditions. 

As can be seen in Figs. |3]j^ bipartite nonlocality and entanglement behave very differently 
in EET in the FMO complex, when the HEOM method is used to model the dynamics. For 
the FRET initial conditions, nonlocality is not found for any pairs of chromophores, while 
entanglement is nonvanishing. Furthermore, for the localized initial conditions, entanglement 
dynamics is similar to the chromophores population dynamics shown in Fig. i.e., oscillating 
and exponentially decaying, while nonlocality does not show any oscillating features at all. 
Instead it drops to zero on a very short time scale, i.e., nonlocality may have a finite life¬ 
time despite the fact that entanglement undergoes an exponential decay. This nonlocality 
sudden death is a phenomenon analogous to entanglement sudden death, i.e., the occurrence 
of finite-time entanglement in coherent systems, which has been predicted^^l and observed 
in quantum optical systems®'. Similarly, nonlocality sudden death has very recently been 
demonstrated in a two-photon experiment® 

CONCLUSIONS 

This work has contributed to the investigation of EET in the FMO complex in considering 
bipartite quantum nonlocality between different chromophores in addition to the bipartite 
quantum entanglement considered in previous studies'^. While entanglement can alterna¬ 
tively be interpreted as measuring the coherences (in the sense of li norm of coherence) in 
the FMO complex, we have shown that nonlocality cannot be given this interpretation, but 
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are given by the diagonal elements pn,..., P 77 of the density operator in Eq. (12). 


should instead be regarded a proper quantiher of quantum correlations between chromophore 
pairs. The numerical simulations show that nonlocality only exists for localized initial con¬ 
ditions. However, it is only observed for those chromophore pairs where one of them receives 
the initial excitation of the system and it drops to zero on a very short time scale (less 
than 100 fs). It should be noted that the behavior of nonlocality found in our simulations 
is related to the fact that the FMO complex can typically exhibit only one excitation at 
a time; by artihcially including more simultaneous excitations would potentially induce a 
much more complicated nonlocality pattern accompanying the EET. The restriction to the 
single-excitation subspace would likely prevent nonlocal multipartite correlations^ among 
more than two chromophores, just as the absence of genuine multipartite entanglement for 
this class of states, well-known from the work by Coffmann et al^. 

The fact that no nonlocality is observed when more realistic initial conditions (FRET 
from the antenna molecules to the FMO) are used indicates that nonlocality is of no im¬ 
portance when considering EET in the FMO in its natural habitat. The entanglement still 
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present, represents local correlations in the sense that they can be described using a theory 
incorporating local realism^. Hence, the correlations between pairs of chroniophores can be 
either quantum or classical, i.e., it cannot be ruled out that EET in the FMO complex just 
as well could be explained from an underlying local realistic framework. Whether this is a 
result of the model used in this study or actually is according to the laws of nature, remains 
an open question. In relation to this latter remark, we further note that to experimentally 
examine the behavior of nonlocality in the FMO complex would be very challenging because 
it would require the ability to perform local measurements in at least two different bases, 
which would be hard to achieve due to the short time scale and the small distances between 
the chromophores. 

Finally, the occurence of nonlocality sudden death found in our simulations with local¬ 
ized initial conditions is another feature that seems to indicate that persistent quantum 
nonlocality is probably rare in biological systems. 
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